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Abstract 

We discuss fluctuation relations in simple cases of non-equilibrium Langevin dynamics. In particular, 
we show that close to non-equilibrium steady states with non-vanishing probability currents some of 
these relations reduce to a modified version of the fluctuation-dissipation theorem. The latter may be 
interpreted as the equilibrium-like relation in the reference frame moving with the mean local velocity 
determined by the probability current. 



1 Introduction 

In statistical mechanics, the fluctuation-dissipation theorem (FDT) provides a simple relation in an equi- 
librium state between the response of the fixed-time averages to small time-dependent perturbations of 
the Hamiltonian and the dynamical correlations 34, [35]. Let O a (x) for a — 1, . . . , A be a collection 
of (classical) observables. With the shorthand notation for the single-time functions O a (xt) of the 
dynamical process x t , the response function and the 2-time correlation function in a steady state are, 
respectively, 

K a "(t-s) = -L-\ {O a t ) and C ab (t - s) = (O a t O b ) , (1.1) 
oris I h=o h 

where ( — ) h denotes the dynamical expectation obtained from the steady state by replacing the time- 
independent Hamiltonian H(x) by a slightly perturbed time-dependent one H(x) — h t O b (x). The FDT 
asserts that when the unperturbed state is the equilibrium at inverse temperature /3 then 

p-^TL^it-s) = d B C ab (t-s). (1.2) 

Such a direct relation between the response and correlation functions is violated in systems out of equi- 
librium and a lot of interest in the research on non-equilibrium statistical mechanics was devoted to such 
violations. In particular, they were studied intensively for glassy systems 13, 10, 7 , for colloidal suspen- 
sions [3:, 19 , for granular matter [2] [6], and for biophysical systems [391 127] . In recent years, it has been 
realized that the FDT, as well as the Green-Kubo relation, another linear response law of the equilib- 
rium regime, are special cases of more general fluctuation relations that hold also far from equilibrium. 
Such relations pertain either to non-stationary transient situations [161 128| or to stationary regimes [20] , 
In particular, the so called Jarzynski equality 28 for the dynamics with a time dependent Hamiltonian 
reduces to the FDT for tiny time variations [8]. 

In the present paper, we revisit the violations of the FDT in simple examples of non-equilibrium 
steady states (NESS) for systems with few degrees of freedom evolving according to the Langevin equation 
possibly including non-conservative forces, see e.g. [421 1141 |4"T1 1251 l23l 144] . For such systems, we show a 
modified fluctuation-dissipation theorem (MFDT) that may be written in the form 

/3-X 6 (*.*) = S.Cf(t,a) (1.3) 

similar to the equilibrium relation, somewhat in the spirit of ref. [44], Above, TZ ab {t,s) and C" 6 (t, s) 
denote the response function and the dynamical correlation function in the Lagrangian frame moving with 
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the mean local velocity v {x) of the NESS. The Lagrangian-frame functions are obtained by replacing in 
the definitions (| 1 . 1 f) the time-independent observables O a (x) by the time-dependent ones O a (t, x) that 
evolve according to the advection equation 

d t O a (t,x) + v {x)-VO a {t,x) = 0, (1.4) 

i.e. are frozen in the Lagrangian frame. In the equilibrium, the mean local velocity v (x) vanishes and 
the MFDT (0 becomes the FDT (fL2)l . 

The other goal of the present work is to explain how the MFDT (|1.3Jl may be obtained from more 
general fluctuation relations by restricting them to the regime close to NESS, similarly as for the case of 
the equilibrium FDT. Before doing that, we recall different fluctuation relations holding arbitrarily far 
from stationarity and equilibrium in Langevin systems and their perturbations. Our discussion follows 
with minor modifications the recent exposition [8]. 

The general results presented in the paper apply, in particular, to two types of one-dimensional systems 
with NESS. The first type describes an exploding Langevin dynamics on the line with a non-Gibbsian 
invariant measure. This case arises, for example, when one studies the tangent process for particles with 
inertia moving in one-dimensional Kraichnan's random velocities |31l 118] . Such velocities vt(x) form a 
Gaussian ensemble with mean zero and covariance 

(v t (x)v s (y))= 5(t-s) D(x-y). (1.5) 

The evolution of the inertial particles is described by the stochastic differential equation (SDE) [4] 

x = u, u = -(— it + v t (x)) , (1.6) 

where r is the Stokes time measuring the time-delay of the particles relative to the flow motion. The 
separation between two infinitesimally close trajectories satisfies then the equations 

^■Sx = 5u, ^8u = i( - 8u + 5x d x v t (x)) , (1.7) 
where one may replace ^d x Vt(x) on the right hand side by a white noise (t with the covariance 

(CtO) = 2r 1 s(t~ s ) (i.8) 

for /T 1 = -^ndlD(Q). For the ratio X = |2, one obtains the SDE 

X = -X 2 - \x + C = -d x H(X) + C (1.9) 

which has the form of a one-dimensional overdamped Langevin equation for the Hamiltonian H(X) = 
\X' A + ±X 2 . The process X t solving Eq. (|1.9Jl escapes in finite time to -co but has a realization with 
trajectories that reappear immediately from +oo, see Appendix A. This corresponds to the solutions 
for (8x, Su) where Sx passes through zero with a non- vanishing speed, i.e. to the crossing of close par- 
ticle trajectories with faster particles overcoming slower ones (allowed in this model of a dilute particle 
suspension with no pressure and no back-reaction on the flow ITT]). The Gibbs density e~' 5£r is not 
normalizable here. The resurrecting process has, however, a non-Gibbsian invariant probability measure 
with constant probability flux. The top Lyapunov exponent for the random dynamical system (|1.6[) is 
obtained as the mean value in this measure of X (which is the temporal logarithmic derivative of \5x\) 
|46| . The above is a variation of a much older story [221 138] of the one- dimensional Anderson localization 
in the stationary Schrodinger equation 

-£r1> + V1> = E1> (1.10) 

with a 5-correlated potential V(x). For Y = (-j^ip)/ip, one obtains the equation 

±Y = -Y 2 - E + V (1.11) 

that may be viewed as a stochastic evolution equation if x is interpreted as time. The invariant measure 
with constant flux for such an SDE was already used in [22] , The substitution Y = X + ^p, E = —j?pz, 
V — C, turns Eq. p. lip to Eq. (|1.9p (provided that one replaces x by t). 

The second particular type of systems with NESS covered by our discussion is obtained by adding a 
non- conservative force to the Langevin dynamics. More specifically, we shall consider a particle which 
moves on a circle according to the SDE 

x = -d x H(x) + G(x) + C, (1.12) 
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where, as before, £ t is the white noise with covariance {L8}. The above dynamics pertains again to 
the overdamped regime where it is the particle velocity rather than the particle acceleration that is 
proportional to the force. The angular coordinate x will be taken modulo 2n. We shall assume that 
H(x + 2ir) — H(x) and G(x + 2tv) = G(x) but G{x)dx / so that the force G is not a gradient 
and it drives the system out of equilibrium. Eq. (|1.12|) was used, for example, to describe the motion of 
a colloidal particle in an optical trap [33] ■ It was discussed recently in 40 in a context similar to the one 
of this work. 

The present paper is organized as follows. Sec. 2, returns to the discussion of stationary Langevin 
diffusion processes, presenting more details on the one-dimensional systems with explicit non-Gibbsian 
invariant measures 8, 40 . For such systems, we examine in Sec. 3 the simplest fluctuation-response 
relation that describes the change of the invariant measure under a small time-independent variation of 
the Hamiltonian. In Sec. 4, we prove the MFDT (|1.3[1 that holds around NESS of the Langevin-type 
dynamics, in particular, in the one-dimensional cases with explicit invariant measures. Sec. 5 is devoted 
to a brief presentation of general fluctuation relations for SDE's [321 1371 24 33 : 8 . These are specified 
for the Langevin systems under consideration in Sec. 6. In particular, we describe the Crooks detailed 
fluctuation relation [12] and the Hatano-Sasa ,24] version of the Jarzynski equality [291 30 , both holding 
arbitrarily far from stationarity. In Sec. 7, we return to the MFDT, showing that it may be viewed as a 
limiting case around the stationary situation of the Crooks transient fluctuation relation or, in a special 
case, of the Jarzynski-Hatano-Sasa equality. Finally, after brief Conclusions, we collect in Appendix A 
some facts about the one-dimensional processes with explosion and illustrate in Appendix B the MFDT 
by an explicit calculation for the Langevin particle driven by a constant force along a circle. 

Acknowledgements. The work of G.F. was supported in part by the US National Science Foundation 
under Grant No. PHY05-51164 and by the Minerva Foundation. 



2 NESS in Langevin processes 

The general stationary dynamics that we consider is described by the Langevin equation in d-dimensions 
with an external force: 

x % = -T ij d j H{x) + IL ii d j H{x) + G'(x)+ (2.1) 

where T is a constant non-negative matrix and LI an antisymmetric one, H is the Hamiltonian, G the 
external force and the white noise £t has the covariance 

(CiCi) = 2p- 1 T i 'S(t-s). (2.2) 

The deterministic force — FX7H decreases the energy, driving the solution towards the minimum of H, if 
it exists, whereas the noise mimics the effect of a thermal bath. We have added the Hamiltonian force 
nV-ff that preserves the energy in order to cover systems governed by Langevin-Kramers equations [32] 
or Fermi-Pasta-Ulam chains [15] • The generator L of the process x t satisfying the SDE (|2.1[) is defined 
by the relation 

a,</(it)> = ((£/)(*«)>. (2.3) 

It is a second order differential operator: 

l = ((-r + n)(Vif) + G) ■ v + /r'v-rv (2.4) 

in the vector notation, with the formal adjoint 

l 1 " = - v ■ ((— l + n)(v#) + g) + /r x v ■ rv . (2.5) 

The transition probabilities Pt(x,dy) of the process satisfy the evolution equation 

d t Pt(x,dy) = L x P t (x,dy), (2.6) 

with the subscript in L x indicating that L acts on the variable x. The dynamics of the mean instan- 
taneous density g t of the process xt is generated by the adjoint operator L* and takes the form of the 
continuity equation 

d t Qt = IS Qt = -V -it (2.7) 

with the currenlQ 

j t = ( - r(Vtf) + n(vff) + g - /r 1 (r - n) v) et . (2.8) 

For convenience, we have included into the current the term f3~ 1 HSJ q. 
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Following [24], let us introduce the mean local velocity vt = Qt^jt- With the use of the velocity field, 
the above continuity equation may be rewritten in the hydrodynamical form as the advection equation 
for the density gt , 

(d t + V-u t )et = 0, (2.9) 

stating that pt is annihilated by the convective derivative or, in other words, that it evolves as the density 
of Lagrangian particles whose trajectories obey the ordinary differential equation 

x = u t {x). (2.10) 



For an invariant density, the corresponding current is conserved: V • j = 0. If for the density g the 
current j itself vanishes then one says that the dynamics (|2,2[l satisfies the detailed balance relative to g. 
The detailed balance holds relative to the Gibbs density e~ l3H if G = (this was assured by the addition 
of the term /3~ 1 Tl\7g to the current). When G / 0, the invariant density is not known explicitly, in 
general, even if it exists. There are, however, special cases of processes satisfying the SDE (|2.2[) where 
one may obtain an analytic formula for a non-Gibbsian invariant density. 



2.1 NESS for resurrecting processes 

The first example, that we shall call type 1 below, is obtained for the Langevin equation on the line. 
In this case, any force is a gradient so that, upon setting for simplicity F = 1, the dynamical equation 
becomes 

x = -{d x H)(x) + C (2.11) 

with the covariance of the white noise given by Eq. (II. 80 . The detailed balance holds here relative to 
the Gibbs density e~^ H since the corresponding current vanishes. Such a density may, however, be 
not normalizable, hence not leading to an invariant probability measure. Let us look closer at various 
possibilities by considering the case of a polynomial Hamiltonian with the highest degree term equal to 
ax k [8]. For k = or k = 1, the solution xt of Eq. (|2.11|l is, up to a linear change of variables, a 
Brownian motion, which does not have an invariant probability measure. For even k > 2 and a > 0, the 
Gibbs measure fx(dx) = Z~ 1 e~ l3H ^dx provides the unique invariant probability measure of the process 
xt- If a < 0, however, then the Gibbs density is not normalizable. In this case the process escapes with 
probability one to ±oo in a finite time and it has no invariant probability measure. For odd k > 3, 
the Gibbs density e~P H ^ is also not normalizable. Changing eventually x to —x, we may assume 
that a > 0. In this case, the process x t escapes in finite time to — oo but it has a realization with 
the trajectories that reappear immediately from +oo, as discussed in Appendix A. Such a resurrecting 
process has a unique invariant probability measure 

X 

fi(dx) = ~(J e 0B(s) dyj e- 0Hix) dx = g H (x)dx, (2.12) 

— oo 

where Z is the (positive) normalization constant. The density g H of fi behaves as (Za[3kx k ~ 1 )~ 1 when 
x — > ±oo, see the estimate (|A.17[l in Appendix A. It corresponds to the constant current j = — (f3Z)~ l 
with the flux towards negative x. The situation provides one of the simplest examples of NESS. In 
particular, the inertial particle in the one-dimensional Kraichnan flow and the Anderson localization in 
the one-dimensional 5-correlated potential lead naturally to the resurrecting processes corresponding to 
k = 3, as discussed in the introduction. 



2.2 NESS for forced diffusions on circle 

The second model with an explicit analytic expression for the invariant non-Gibbsian measure, that we 
shall call type 2 below, is the perturbed one-dimensional Langevin equation (I1.12|l on the unit circle. 
Now, the unique invariant probability measure is given by the formula [421 1251 14U] 

^dx) = e mx ' v) dy) e- 011 ^ dx = g H (x)dx (2.13) 



for < x < 2tt, with 



X X A-K 

U(x, y) = H{y) + Q(x - y) J G(z) dz + 6{y - x) ( j G(z) dz + J G(z) dz^ . 



(2.14) 
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Also here the density g H of the measure n corresponds to a constant probability current 

7 = ^(e - -l). (2.15) 

In the following, we shall see how the presence of the probability flux in NESS deforms the usual fluctuation 
relations. 

3 Modified fluctuation-response relation 

As a warmup, let us see what is the form taken by the most elementary fluctuation-response relation 36 in 
the one-dimensional systems with NESS that we discussed above. The setup of the fluctuation-response 
relation is as follows. One prepares the system in the far past in the invariant state with probability 
density g H that is given by Eqs. (|2.12|l and (|2.13|l for the type 1 and type 2 systems, respectively. At 
t = 0, the Hamiltonian H is perturbed by a small time-independent potential V (vanishing sufficiently 
fast when x — * ±oo for type 1), leading to the change H ^ H' = H + V . The systems evolves then 
and converges toward the new steady state with the probability density g , . The fluctuation-response 
relation compares the initial and the final averages of an observable Ot = 0(xt) '■ 



(O ) = J 0(x)g H (x)dx and (<?«,) = J 0(x) g H , (x) dx . (3.1) 

By a straightforward differentiation of the explicit formulae for the invariant densities, one obtains the 
identity 

<<U = (O > -P[(O (V - V )) - (O )(V - %)] , (3.2) 
up to terms of the second order in V, where 

j V(y)e^ H ^dy f V(y) e 0u ^dy 

V(x) = ^ and V(x) = ^ (3.3) 
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for the systems of type 1 and type 2, respectively. Eq. (|3.2p deforms the usual fluctuation-response relation 
around the Gibbs state by replacing V by (V — V) with the averaged potential V dependent on 
the initial Hamiltonian H . 

4 Modified fluctuation-dissipation theorem 

Coming back to the general case, let us consider a system prepared at negative times in the steady state of 
the stationary Langevin dynamics (|2.1|l . This forces the time zero value x of the corresponding process 
to be distributed according to the invariant probability measure (x (dx) — g (x) dx. At t — 0, one switch 
on a non-stationary perturbation taking the Hamiltonian for the positive times to be equal to 

H t (x) = H{x) - J2 h «,tO a (x) , (4.1) 

a 

where h a ,t carry the time dependence and functions O a (x) (the "observables" ) are supposed, for sim- 
plicity, to vanish sufficiently fast when \x\ — » oo. We denote by {j^ the corresponding expectation, 

with (jF) referring to the non-perturbed case. The expression 

{TRt) =7^-1 {T). (4.2) 

defines the response correlations. To shorten further the notations, let us set 

{O a t R h s \ = ll ab {t-s), {O a t O b s ) a = C ab (t-s), e(t-s)(O?B b s ) = B ab (t-s) (4.3) 

for O t = 0(x t ) and the induced observable 

B" = q~\ ■ VO" (4.4) 

with j standing for the current H2.8[) corresponding to the invariant density g . Note that for the 
one-dimensional NESS with constant probability current j , the observable B b — j gg 1 d x O b has the 
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probability flux as an explicit factor. Remark that, by causality, the response function 1Z ab (t — s) vanishes 
for s > t. We shall show the following modified version of the fluctuation-dissipation theorem (MFDT) 
holding for t > s: 

/T 1 Tl ab {t - s) = d s C ab (t - s) - B ab (t - s) . (4.5) 



The second term on the right hand side of Eq. (14.5[) is new as compared to the FDT around the equilibrium 
steady state. Indeed, in the equilibrium case, the external force G = and g = Z~ 1 e~ ,3H is the 
normalized Gibbs factor with j = so that B b = and the MFDT (|4.5[) reduces to the standard 
equilibrium form (|1.2[l , 

4.1 Lagrangian-frame interpretation 

What the formula (|4.5|l for the response function means is better understood by rewriting it with the 
explicit form of the right hand side as 

p-'K^it-s) = Jdx O b (x)g (x) d s P t - s (x,dy) O a (y) -Jdx (j -S7O b )(x) P t - S {x,dy) O a {y) (4.6) 

r dx O b (x) (d a + V-u a (x)) f g (x)P t - s (x,dy) O a (y), 



where Pt(x,dy) denotes the stationary transition probabilities of the unperturbed process and we have 
integrated by parts to obtain the second equality, setting v = Qg l j - Note that the time derivative d s 
of the equilibrium relation has been replaced by the convective derivative d s + V ■ v Q (x) which acts on the 
first component of the joint probability density function of the time s and time t values of the stationary 
process xt- This suggests that the MFDT (|4.5p should take the equilibrium form in the Lagrangian frame 
moving with the stationary mean local velocity v Q (x). 

To render this interpretation more transparent, let us replace the time-independent observables O a (x) 
by the time dependent ones O a (t,x) evolving according to the advection equation (|1.4|l . We shall define 
the Lagrangian-frame response function and correlations function by 

^(*' S ) = ^-| <°?(*)> h ' C a L b (t,s) = (OUt)O b (s)) (4.7) 
where ( — ) denotes now the expectation referring to the Hamiltonian Ht(x) = H(x) — ^ h t , a O a (t, x). 



Ot(t) = 0(t,x t ) with the double time-dependence. Writing the MFDT (|4.5[) with the explicit right hand 
side given by Eq. (14. 6[) for the observables O a (t,x), one casts this relation into the form 



ff-'K" M) = / dx O b (s,x) Qo (x) d a Pt-s(x,dy) O a {t,y) 

dx g (x) (v {x) ■ VO b (s,x)) Pt- S (x,dy) O a (t,y) 
= 8 S I dxO b (s,x)g (x) Pt- S {x,dy) O a {t,y) = d s C?(t,s), (4.8) 



where the last equality follows from the advection equation (|1.4[) . This proves the identity (|1.3[) announced 
in the introduction. We should stress that, in spite of the similarity between that relation and the 
equilibrium FDT (|1.2p . in general it is not true that the dynamical process x t viewed in the Lagrangian 
frame of the velocity field uq is governed by an equilibrium Langevin equation, although this is what 
happens in the simple example considered in Appendix B. 

For the Langevin process on the circle, a fluctuation-dissipation relation for velocities similar to (14.51) 
was discussed in [44], see Eq. (11) there, with the interpretation similar in spirit, but not in form, to the 
above one, see the subsequent discussion there. One of the consequences of the fluctuation-dissipation 
relation of 44 linking the effective diffusivity and mobility was checked experimentally in [5], see also 

M- 

It is sometimes more interesting, especially for applications, to re-express the fluctuation-dissipation 
relations in terms of the frequency space quantities. Let 



oo oo 



C ab (io) = J e l " ( *- s) C ab (t- s)dt = J e luj(t - a) C ab (t- s)dt + J e" l " (t - s) C ba {t - s) dt , (4.9) 

— OO s s 
oo oo 

H ab {iu) = I e lt " ( '- s) TZ ab (t -s)dt, B ab {uj) = I e llj(t - s) B ab {t - s) dt . (4.10) 
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Note that TZ ab (u>) measures the response to the time-dependent potential of frequency uj. The MFDT 
()4.5|) is equivalent to the relation 

iud a \u>) = /3- 1 (H ab (uj) -1l ba {-Lu)) + B ab {uj) - B ba {-uj) (4.11) 

in the frequency space. 

In general, assuming that the transition probabilities Pt{x, dy) converge at long times to the invariant 
measure, all three terms of Eq. (|4.5jl tend to zero when (t — s) — > oo. Mimicking the idea employed with 
success for disordered systems |13l 1141 [7], their relative proportions, or the relative proportions of the 
corresponding terms in Eq. (|4.11|l , could be used to define dynamical temperatures that would, in general, 
depend also on the observables involved. We discuss those proportions in a simple case of the Langevin 
dynamics on a circle with a constant force in Appendix B. 

We propose three derivations of the result (|4.5[l : the first one direct, that we shall present now, and 
the next two ones from the general fluctuation relations that will be discussed in the subsequent sections. 



4.2 Direct derivation 

The beginning of the argument is quite standard, see e.g. [T] or Sec. 2.3.2 of |41] , By the definition of 
the response correlations, 

K ab (t-s) = [o a (y)g t (y)dy, (4.12) 

where gt is the density obtained by the perturbed dynamical evolution (|2.7|) from g . Using the explicit 
form (|2,8[l of the current, one obtains by the first order perturbation the relation 

t 

Qt(y)dy = g a (y)dy-J2 f h b , a ds j dx [V • ((r - n)(V0 6 ) g )] (x) P t _ a (x, dy) + 0(h 2 ) . (4.13) 
b o 

Consequently, for t > s, 

1Z ab {t-s) = -jdx[V-{(r~U)(VO b )g )](x)P t ^(x,dy)O a (y). (4.14) 
Now, a straightforward although somewhat tedious algebra shows that 

/r'v • ((r - n)(vo 6 ) Bo ) = L\o b Ba ) + v ■ (o b j a ) = L\o b Ba ) + j a ■ vo b , (415) 

where the adjoint generator is given by Eq. (|2,5[l and in the last equality we have used the conservation 
of the current j . Substituting this identity to Eq. (|4.14|l and integrating the term with U by parts, we 
obtain the relation 

j3- x Ti ab {t-s) = - J dx[(O b g„)(x)L x + (j -VO b )(x)]P t . a (x,dy)O a (y) (4.16) 
which, together with Eq. lf2~6) , implies the MFDT (P). 



5 General fluctuation relations 

In [S], two of us discussed arbitrary diffusion processes in d dimensions defined by the Stratonovich SDE 

x — u t (x) + v t (x) , (5.1) 

where Ut(x) is a time-dependent deterministic vector field (a drift) and v t (x) is a Gaussian random 
vector field with mean zero and covariance 

(vl(x)vl(y)) = 6(t-8)D?(x,y). (5.2) 

The Langevin equation (|2.1[) provides a special example of such an SDE. For the processes solving 
Eq. (|5.1[) . we showed, combining the Girsanov and the Feynman-Kac formulae, a detailed fluctuation 
relation (DFR) 

» (dx) P T (x-dy,dW) e~ w = ii' {dy*) P^(y*;dx* ,d(-W)) , (5.3) 

where 

1. ^ (dx) = g (x)dx is the initial distribution of the original (forward) process, 
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2. y,'(dx) — g' (x)dx is the initial distribution of the backward process obtained from the forward 
process by applying a time inversion (see below) , 

3. P T (x; dy, dW) is the joint probability distribution of the time T position x T of the forward process 
starting at time zero at x and of a functional W T of the same process on the interval [0, T] 
(described later), 

4. P T (x;dy,dW) is the similar joint probability distribution for the backward process. 

The key behind the DFR is the action of a time inversion on the forward system. First, such an inversion 
acts on time and space by an involution 

(t,x) i — ► (T-t,x*). (5.4) 

The above involution may be extended to the action x x on trajectories by the formula xt = x^-t 
and, further, to the action on functionals of trajectories T h- ► T by setting J-[x] = F[x\. Second, to 
recover a variety of fluctuation relations discussed in the literature .32 , 37, 11, 12, 28, 45, 9], we allow for 
a non-trivial behavior of the drift ut(x) under the time inversion, dividing it into two parts, u = u+ +u_, 
with u+ transforming as a vector field under the space-time involution ()5.4[) and u~ as a pseudo-vector 
field, i.e. defining 

u T-t,±{x*) — ±{d k x* l )(x) u t> ±{x) , u — u' + + u'_ . (5.5) 

The random field v t (x) may be transformed with either rule. By definition, the backward process satisfies 
then the Stratonovich SDE 

x = u' t (x) + v' t (x) (5.6) 

and is, in general, different from the naive time inversion xt of the forward process. The functional W T 
involved in the DFR depends explicitly on the densities g a and g T , where the latter is defined by the 
relation fi' (dx*) = q t (x) dx : 

T 

W T = -A T \ng + jjtdt (5.7) 

o 

with the notation A T In g = In g T (x T ) — In g (x ) . In the above formula, 

Jt = 2ut,+ (xt)-dt 1 {x t )(xt-ut,-(xt)) -(V •«*,_)(»«), (5.8) 

where d\ 3 (x) = D\ 3 {x,x) and u\ t +{x) = u\ + (x) — \d j D\° (x, y)\ v=x . The time integral in Eq. (15. 7p 
should be taken in the Stratonovich sense. The functional W' T for the backward process is defined in 
the same way, setting fj, {dx*) = g' T (x)dx. One has the relation 

W T = -W T , (5.9) 

where the tilde denotes the involution of trajectory functionals introduced before. 

The quantity J t has the interpretation of the rate of entropy production in the environment modeled 
by the thermal noise. When the density g T coincides with the density obtained from £ by the dynamical 
evolution (|2.7p . where now the current 

jt = (ut - ^d t \7)g t with u = u+ + it- , (5.10) 

then the first contribution — A T In £> to W T may be interpreted as the change in the instantaneous 
entropy of the process. In this case, the functional W T becomes equal to the overall entropy production. 
Keep in mind that this is a fluctuating quantity which, in general, may take both positive and negative 
values. 

The DFR (|5.3[) holds even if the measures jj, and fj,' B are not normalized, or even not normalizable. 
For normalized initial measures, we denote by 

(F) = J T[x\M[dx] and (F) r = J F[x] M'[dx] (5.11) 

the averages over the realizations of the forward and the backward process xt, < t < T, with x 
distributed according to the probability measure fj, and //, respectively. A^[da;] and _M'[d:r] stand 
for the corresponding measures over the space of trajectories. One of the immediate consequences of the 
DFR JO) is the identity 

<e- w ^) = 1 (5.12) 

obtained by the integration of the both sides. This is a generalization of the celebrated Jarzynski equality 
[29 1 I28|. The relation (|5.12[) implies the inequality (W T ) > that has the form of the second law 
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of thermodynamics stating the positivity of the average entropy production. The Jarzynski equality, 
however, provides also information about an exponential suppression of the events with negative entropy 
production in non-equilibrium systems for which (W T ) > 0. 

With a little more work [8| involving a multiple superposition of the FDR (|5.3[) . the latter may be 
cast into the Crooks form [T2] 

(Tc- w t) = (;F}'. (5.13) 
In terms of the trajectory measures .M[da;] and A1'[da;], this becomes the identity 

M'[dx] = e- w T^M[dx], (5.14) 

where M'[dx\ = A^'[daT|. Eq. (|5.14|l permits to interpret the expectation of W T as the relative entropy 
of the trajectory measures: 

<W T ) = fhx M[dx] = S(M\M'), (5.15) 

J M'[dx] 

in line with the above entropic interpretation of the functional W T . 



6 Fluctuation relations in Langevin dynamics 

Let us specify the DFR (|5.3[) to the case of Langevin dynamics (|2.2|l with, possibly, time-dependent 
Hamiltonian H and external force G. A canonical choice of the time inversion for such a system takes 

u+ = -TVH, u- = ILVH + G (6.1) 

and a linear involution x* = Rx such that RTR 1 — F and RUR 1 = —II 8 . For example, for 
the Langevin-Kramers dynamics in the phase-space, the usual R changes the sign of momenta. The 
backward dynamics has now the same form as the forward one, with the time-reversed Hamiltonian 
H' t (x) = Hr-t(Rx) and the time-reversed external force G' t (x) = —RGr-t(Rx). In this case, we may 
use the Gibbs densities Qt(x) = Z^ 1 e~ /3Ht ^ = g' T _ t (Rx) at the initial and final times, with Zt standing 
for the partition function J e _l3Ht ^ x 'dx if the integral is finite and Zt = 1 otherwise. A straightforward 
calculation gives: 

T 

W T = ln(Z T /Z a ) + J {(5dtH t + f3G t ■ VH t - V ■ G t )(x t ) dt . (6.2) 
o 

For normalizable Gibbs factors this is often called the "dissipative work". W^, is given by the same 
expression with Ht and Gt replaced by H[ and G' t . 

Another useful choice of the time inversion is based on the eventual knowledge of the densities gt 
corresponding to the conserved currents with V • jt = 0. Note that such densities would be left invariant 
by the evolution (|2,7[1 if the time-dependence of the Hamiltonian and of the external force were frozen to 
the instantaneous values H t and Gt- One takes 

u+ =/3 -1 rvine, u- = -rv{H + p^ing) + nv# + g. (6.3) 

With the linear involution x* — Rx as above, the backward process has 

u+ = /ravine', u'_ = rv(H' + /r 1 in £ /) + nv#' + g' , (6.4) 

where g't(x) — gr-t(Rx) and H' and G' are as before. The current corresponding to the density g' t 
satisfies 

j' t (x) = -Rj T _ t (Rx). (6.5) 

It is also conserved. Such a time inversion (for 7? = 1) was considered in [23] and, more explicitly, in |S]. 
In 0, it was called the current reversal. The DFR (|5.3[) holds now for 

T 

W T = - J (dtlng t )(xt)dt (6.6) 
o 

and given by the same expression with gt replaced by g' t . The Jarzynski equality (|5.12|l for this case 
(assuming that g t are normalized) was first proven by Hatano and Sasa in [24 . Note that if G — then 
the current corresponding to the densities g t = Z t _1 e~ /3Ht is conserved and with this choice of gt the 
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two time inversions coincide. In particular, for G = and the time-independent H t = H , the functional 
W T identically vanishes and the DFR reduces to the equality 

e~ pB( - x) dx P T (x,dy) = e - 0H{y) dy P T {y* , dx") (6.7) 

which is a more global version of the detailed balance relation. On the right hand side, one may replace 
P T by P T for II = and x* = x because in that case, the forward and backward processes have the 
same distribution.. 



6.1 Fluctuation relations for resurrecting processes 

Let us consider the process solving the SDE (|2.11[) , with H t (x) = ax k + o(\x\ k ) at large \x\ with odd 
k > 3 and a > 0. We admit a mild time-dependence of Ht(x) disappearing when x — > ±oo. The 
corresponding process still has a resurrecting version, as in the stationary case described in Appendix A. 
Let us discuss first the canonical time inversion with the trivial involution x* — x leading to the backward 
process of the same type with Ht(x) replaced by H' t (x) — HT-t(x). The definition of the functional 
W T for the resurrecting process requires a little care in order to account for the contributions from the 
jumps from — oo to +oo. This may be done by compactifying the line to a circle writing x = cot 9 for 
6 modulo 7r. One has 

XT T 

J J t dt = -P J (d x H t )(x t ) x t dt = j {(3ka cot*" 1 t + . . . ) sin" 2 6 t 6 t dt (6.8) 



and the integral diverges to +oo whenever 6 t passes from the negative to the positive values, i.e. 
whenever x t jumps from — oo to +oo. Upon taking g t = e~ 0Ht , we infer that 

T 

.J (dtH t )(xt) dt if x t has no rebirths for < t < T , 

Vv T = \ o (6.9) 

+oo otherwise . 

and similarly for W' . As a result, the contributions of the rebirths to the DFR (|5.3|) trivially decouple 
reducing the latter to the identity 

e" Ho{x) dx P°{x;dy,dW) e' w = e~ H T^ dy P^° (x; dy,d(-W)) (6.10) 

between the joint distributions of the endpoints and of the functional W T — (dtHt)(x t ) (usually called 
the "Jarzynski work"), or of its counterpart W~, in the processes without rebirths. In the stationary 
case, Eq. (|6 . 10 p reduces to the detailed balance relation 

e - 8H(x) dx P?(x,dy) = e- 0HM dy P?(y,dx). (6.11) 

which assures upon integration over x that the process without rebirths preserves the infinite measure 
e- 0H ^dx. 

On the other hand, one could use for the same SDE with the resurrecting solution the current reversal 
based on the splitting 

u t ,+ = P~ 1 d x \riQ Ht , ut,- = -d x {H t + In g H J (6-12) 

of the drift —d x Ht, with the density g H given by Eq. (|2.12|l . The use of the involution x* = —x leads 
then to the backward process solving the SDE (|2.11[) with the Hamiltonian Ht(x) replaced by 

—x 

H' t (x) = H T _ t (-x) - 2 / 3" 1 ln ( f e f3H T-^dy^j . (6.13) 

— oo 

From the estimate (| A. 17[) in Appendix A, one infers that H' t (x) — ax k + o(x k ^ 1 ) for large \x\ (we 
have used the non-trivial spatial involution to keep a positive). Hence Ht is of the same type as the 
Hamiltonian Ht for the forward process. In this case, the functionals W T and are given by Eq. (|6.6p 
with Qt(x) replaced by Q Ht (x) and q Ht ^—x) = q h ,(x), respectively, with no extra contributions from 

the rebirths. In the stationary case, one has W T = = and the DFR (|5.3[) reduces to the modified 
detailed balance relation 

q h (x) dx P T (as, dy) = g H (y) dy P' T (-y, d(-x)) . (6.14) 

The latter links the transition probabilities of the resurrecting forward and backward processes and 
assures upon integration over x or y that those processes preserve the probability measures Q H (x) dx 
and Q H (—x)dx, respectively. 
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6.2 Fluctuation relations for forced diffusions on circle 



Consider the process satisfying the SDE (|1.12|l with periodic Hamiltonian H(x) — H(x + 2ir) and external 
force G(x) = G(x + 2n), both possibly time-dependent. The use of the canonical time inversion with 
the trivial involution x* — x leads to the backward process of the same type with H' t = H' T _ t and 
G' t = — G' T _ t . The functionals W T and Wj> are given here by the formula (|6.2[) . 

On the other hand, the use of the current reversal with the densities g H of Eq. (|2.13p and the trivial 
inversion x* — x leads to the backward process satisfying the same SDE with H t replaced by 

H' t {x) = -H T _ t (x) -2/3- 1 lng HT _ t (6.15) 

and Gt by G' t as above. The functionals W T and are given now by Eq. (|6.6p with gt(x) equal 
to g Ht (x) and g H ^(x) = g H (a;), respectively. They vanish in the stationary case when the DFR (15. 3[) 
reduces again the modified detailed balance relation 

g H (x)dx P T (x,dy) = g H (y)dy P^(y,dx) . (6.16) 

Note that, in general, to obtain the DFR for the probability distributions on the circle, one should sum 
both sides of the relation (|5.3[1 pertaining to the motion on the line, over the shifts of x or y by 2-nn with 
integer n (both summations amount to the same). To get the Jarzynski equality, one has to integrate 
the relation obtained this way over x and y from to 2n. Another simple remark is that for the SDE 
(|1.12[) on the circle, one may always assume that the external force G is constant by changing G(x) to 
G = 27 Jo^ G(x)dx and by subtracting f*(G(y) — G)dy from H(x). Such a change does not affect 
the DFR (|5.3[) obtained by the current reversal that uses only the invariant densities and it modifies in a 
simple way the DFR obtained from the canonical time inversion because in the latter we used the Gibbs 
measures for the initial and final distributions. 



7 Fluctuation relations close to NESS 



As promised, we shall show here that the MFDT (|4.5[1 . proven directly in Sec. 4.1, may be also derived 
by reducing the Crooks DFR for the current reversal, and, in a special case, the Jarzynski-Hatano-Sasa 
equality, to the situations close to NESS. 

7.1 Reduction of Crooks' DFR to MFDT 

Let us consider the DFR (|5.13[1 for the Langevin dynamics (|2.1I) with T = O™ and < t < T, the 
backward dynamics determined by the current reversal, and, for simplicity, the trivial space involution 
x* = x, see Sec. [6] It reads: 

(O a t e~ W T) = (O a T _ t )'. (7.1) 



We shall assume that the time-dependence of the Hamiltonian is given by Eq. (|4.1|l with h a ,o = = 
h a ,T, and that the external force G is time-independent. Let, as above, g (x)dx denote the invariant 
probability measure of the unperturbed process (assumed to exist). By gt, we shall denote now the 
normalized densities whose current j t , given by Eq. (|2.8[) . is conserved, i.e. such that 



tfet - Yl h ^ v ■ (( r - n )(vo a ) gt) = o . (7.2) 

a 

Expanding 



= Go + + ( 7 ' 3 ) 



we obtain from Eq. (|7.2[) the relation 

L^el = v- ((r-n)(vo a ) ft) ) (7.4) 

whose unique solution 

el = {tf)- x [V-({T-Tl){VO a )e )} (7.5) 

is chosen by imposing the orthogonality of gf to the constant mode, required by the normalization of gt. 
For the current reversal, the functional W T is given by Eq. (|6.6|l so that 



H> = - 



1 

J (9tht, a ) {Q^Ql) t dt + 0(h 2 ) 
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T 

= E/ [v ■ ((r - n)(vo a ) a,)] } t dt + o(h 2 ) , (7.6) 

a Q 

where we have integrated ones by parts over t. The application of the operator s £ \h=o for < s < t 
to the both sides of Eq. (|7.1[) gives the identity 

n a \t~s) - (o a t d s { e -\i))- x [v ■ ((r -n)(vo a ) Qa )]} s ) = o. (7.7) 

We have used the fact that, by causality, the right hand side of Eq. (|7.ip does not give the contribution 
because the perturbation is concentrated around time (T — s) > (T — t). From the relation (|7.7[l . it 
follows that 

n ab (t-s) = (o t a 9 s { f?0 - 1 (L t )- 1 [v-((r-n)(vo a ) eo )]} s ) 

= 8. J dx {(L^) -1 [V ■ ((r — II)(VO a ) g )] }(x) P t - S (x, dy) O a (y) 

= - J dx {(Lt)- 1 [V ■ ((r - n)(VO a ) g )]}{x) L x P t - s {x, dy) O a (y) 

= - j dx [V ■ ((r - n)(VO a ) ft,)] (x) P t _ s (x, dy) O a (y) (7.8) 

which is Eq. (|4.14p above. The rest of the proof of the identity (|4.5p goes as before. 



7.2 Jarzynski-Hatano-Sasa equality and MFDT 

The standard FDT around the equilibrium Langevin dynamics (|2.ip without the external force may be 
obtained by expanding the Jarzynski equality (|5.12[) for the Hamiltonian (|4.1|l up to the second order in h, 
see [8]. The Jarzynski equality may then be viewed as an extension of the FDT to the case of Hamiltonians 
with arbitrary time dependence driving the system far from equilibrium. The natural question is whether 
this picture may be generalized to the case of the modified FDT (|4.5|) holding around NESS. We shall 
show here that the answer is a qualified yes. 

Let us expand to the second order in h the Hatano-Sasa version of the Jarzynski equality (|5.12p 
obtained from the Croocks DFR (|7.ip for the current reversal by replacing O a by 1. We shall need to 
know the form of the densities Qt with conserved current, i.e. satisfying Eq. (|7.2p . to the second order in 
h. One has 

Qt — Qo + E^*' a£l i + E ht,ah t ,bQ2 + 0(h 3 ) (7.9) 

a a.b 

with Qia as before, see Eq. (|7.5|L and 

e f = (L+r 1 [v ■ ((r - n)(vo°) q\)} . (7.10) 

Expanding, in turn, the functional W T given by Eq. (|6.6p to the second order, we obtain 

T T 

w T = -Y,j (fcMOrt)**"!! / (WKaKb)) {q^q? - te; 2 eiQ b i) t dt + o(h 3 ), (7.11) 
a a ' b 

and further, 

T T 

(e~ w -) - 1 = E/ (*M((rt)t)*+E / (ft(h*..^))((rt b -5ei- a e?^) t ) * 

a a,6 Q 

T T 

+ iE/ f(dtK«){dsh s , b ) ({Q- 1 eD t {Q^Q\) s ) o dtd a + 0(h 3 ). (7.12) 

The second term on the right hand side integrates to zero because of the stationarity of the unper- 
turbed expectation and the boundary conditions ho, a = = hr, a - Expanding the remaining perturbed 
expectation in the first term on the right hand side, we infer that 

T T 

(e" W ^) - 1 = E/ j{dth t , a )h s , b ^{e- 1 el) t R b \ dtds 

a,b Q Q 
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T T 

+ f (dM (d s h s , b ) (( e ; 1 Qt) t fa 1 Q b 1 ) B ) o dtds + 0(h 3 ) = 0, (7.13) 

a,b o 

where the last equality follows from the (generalized) Jarzynski equality (|5.12|) . The integration by parts 
and the causality permit to conclude that for t > s, 

d t ((0-V) t -R') o = 9 t 9 s ((f? r 7 1 rf) t (f? r 7 1 f?i) s ) o , (7.14) 
or, integrating ones over t, that 

= d a ^(g~ 1 Qi) i (Q~' L Qi) a j (7.15) 

(we used the fact that both sides vanish for t = s). Note that Eq. (|7.15[) stays true if we add to gl any 
multiple of go, so that we may drop the normalization condition J Qi(x) dx — 0. From Eqs. (|7.5|l and 
(|4.15[) . it follows then that we may take 

Q^Qi = P{O a -d a ) = (1A\ (7.16) 

where the dressed observable 

6 a = -Q-\l))- l [j -VO a }. (7.17) 
The identity (|7.15|l becomes now the relation 

p-\A?Rl) = d s (A a t Al) = d s (AZO b s ) - d s J dx (d b g )(x)P t . s (x,dy)A a (y) 

= d s (A a t O b 3 \ + j dx (d b g )(x) L x P t ^(x,dy)A a (y). (7.18) 

Integrating by parts in the last term and using the definition (|7.17[) . we finally obtain the identity 

p-\A1 R b \ = d s (A a t O b ) Q - j dx (j ■ VO b )(x) P t „ s {x,dy) A a (y) 

= d 3 (A?O b s ) - (A?B b a ) Q (7.19) 

which is the MFDT (|4.5|l with the observable O a replaced by A a . It is a consequence of the identity 
H4.5|) but, in general, it does not seem to be equivalent to it, except for the equilibrium case with vanishing 
external force and g = Z~ 1 e~ 0H when j = and A a = O a . 

In the special case of the one-dimensional NESS with constant current j described above, the dressing 
(|7.17[) of the observables coincides with the one given by Eqs. (|3.3|l for the types 1 and 2, respectively. 
This may be easily seen by checking that for the latter, 

L f Q V=-j d x V (7.20) 

with g = g H given by Eqs. (|2.12|l and (|2.13p . 



8 Conclusions 

We have discussed different fluctuation relations for the Langevin dynamics. Those included an extension 
of the fluctuation-dissipation theorem (FDT) (|1.2[) . one of the most important relations of the (close 
to) equilibrium statistical mechanics, to the case of non-equilibrium steady states (NESS) of Langevin 
processes. The modified fluctuation-dissipation theorem (MFDT) ()4.5|) that holds around NESS has a 
new term containing the probability current but in the Lagrangian frame moving with the mean local 
velocity determined by the current, it takes the form (|1.3p similar to that of the equilibrium FDT. We 
also pointed out that, similarly to the equilibrium FDT, the MFDT may be viewed as a limiting case 
of more general fluctuation relations that are valid arbitrarily far from the stationary situation, namely 
of the Crooks detailed fluctuation relation (|5.13p for the backward dynamics with inverted probability 
current or of the Jarzynski-Hatano-Sasa equality (|5.12p . The general discussion was illustrated on two 
examples of one-dimensional systems with explicit non-equilibrium invariant measures. 
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Appendix A 

We return here to the case of the stationary SDE (|2.11|l . The transition probabilities Pt(x,dy) for 
the diffusion process solving this equation are given by the kernels of the exponential of the generator 
L = f3~ 1 d x — (d x H)d x of the process: 

J P t (x,dy)f(y) = (e tL f)(x), (A.l) 
see Eq. (|2,6|l . One has the following relation 

L = -r^^Q^Qe-i 011 (A.2) 
for Q = d x + \(3{d x H), Q 1 = -d x + \(3{d x H). The Fokker-Planck operator 

pr^OHQ = -P^dl + \{d 2 x H) + j(3(d x H) 2 (A.3) 

is a positive self-adjoint Hamiltonian of a super-symmetric quantum mechanics |47] so that the transition 
probabilities may be defined by the relation 

Pt(x,dy) = e^" W e- tp ~ lQ ' Q {x,dy) e "i™ (A.4) 

If the Gibbs density is normalizable with Z = J e~ H ^dx < oo then ip (x) = Z~ 1/2 e~^^ H ^ provides 
the zero-energy groundstate of the Fokker-Planck Hamiltonian. Such a groundstate is supersymmetric: 
Qipo = 0. For the Hamiltonian H(x) = ax k + o(|x| fc ) at large |x| with either even k > 2 and a < or 
odd k > 3, however, the groundstate V °f /3~ 1 Q t Q is not given by eT^I 311 ^ but has positive energy 
E > and breaks the supersymmetry: Qi{> ^ 0. In these cases, the transition probabilities (|A,4|) are 
not normalized for t > with 

1 > / Pt(x,dy) ~ e"^*. (A.5) 

/ £ large 

The defect (l - / P t (x, dy)) gives the probability that the diffusion process xt solving the SDE (|2.11[) 
and starting at time zero at x escapes by time t to ±oo. Writing Pt(x, dy) = Pt{x, y) dy, the probability 
that the escape happens between times s and s + ds may be expressed as 

(-£ J P s (x,dy))ds = (- j ' d y [{p- 1 d y + {d y H){v))P s {x, V )\dy)ds 

= (-(r'dy + idyHUy^Psix^^Jjds (A.6) 

with the y = ±oo terms determining the rates of escape to ±oo, respectively. 

Let us concentrate on the case with H(x) = ax k + o(|a;| fc ) for odd k > 3 and a > 0, denoting 
the transition probabilities of Eq. (|A.4|) by Pt{x,dy). Although they are not given by a closed analytic 
expression, their time integral, equal to the Green kernel of L, is: 

oo min(x,y) 

J dsP°{x,dy) = (~L)~ 1 (x, dy) = /?( J e 8H(z) dzj <T SHbl) dy . (A.7) 

-oo 

From the last formula, we infer that 

oo 

lim fds(p- 1 dy + (d y H)(y))P°(x,dy) = 6 1 ^ 1 (A.8) 

y^±oo J 


so that the process xt escapes here only to — oo (this is already true if we ignore the noise in Eq. (|2. 11 [I ) . 
On the other hand, since the limit of the right hand side of Eq. (|A.7|) when x — » +oo exists, it follows 
that the transition probabilities from x = +oo, 

P?(oo,dy) = lim P?{x,dy) (A.9) 

X — > + oo 

are finite, non-zero measures. One may then define a resurrecting version of the process x t solving the 
SDE (|2.1ip . The trajectories of such a process, after almost surely reaching — oo reappear immediately 
at +oo. The resurrecting process is Markov and its transition probabilities are 

oo 

Pt(x,dy) = ^P^(x,dy), (A. 10) 
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where P^{x,dy) are the transition probabilities with exactly n jumps from — oo to +00. They are 
given by the recursion relation: 

t 

P t n+1 (x,dy) = -Jds(±j P°(x,dy))Pr- s (oo,dy). (A.ll) 


In terms of the Laplace transforms 

00 

P"(x,dy) = J e- tu P t n (x,dy), (A.12) 


the recursion (|A.11|I becomes the equality 

PZ +1 (x,dy) = (l-u J P°{x,dz))P^ (00, dy) (A.13) 

which may be easily solved by iteration: 

PZ(x,dy) = (l-u J P°(x, dz)jp° (00, dy) (A.14) 

for n > 1. Re-summing the geometric progression, one obtains for the Laplace transform of the transition 
probabilities of the resurrecting process the expression 

(l-wj P°(x,dz))P° (00, dy) 

P^{x,dy) = P°(x,dy) + ± ' . (A.15) 

w J PiJ(co, dz) 

Note that P®(x, dy) is analytic in uj for Recj > — E Q but P u (x, dy) has a pole at u> — with the residue 

— oc 

where Z is the normalization constant. This is the invariant probability measure (|2.12|l of the resurrecting 
process. 

Let us finish by estimating the behavior of the density of the invariant measure fJ,(dy) when \y\ — > 00. 
We shall show that 

| e ,(H(«)) rfz= _^ + o( ^ +1) _ (A1?) 

To this end, we rewrite the latter integral as 

00 00 

J e -?my)-«(y-z)) dz = y i-k J e - BWy) -H(y-uv^)) du (A18) 



We take H (y) = ay h + h(y) and assume that h(y) is smooth and that 

Kv)\y\- k , 0, (dyh)(y)\y\- h+1 -» 0. (A.19) 

First note that for z > 0, 

k 

a(y k ~(y~z) k ) = a]T(-l) i+1 C4 y k ~ l z l > Zezy^ 1 + ez k (A.20) 

if e > is small enough. Next, for < z < ||y| and \y\ big enough, 

|fc(y) - h{y - z)\ < z\d y h(y - $z)\ < ezy^ 1 (A.21) 

and for z > ||y|, 

\h(y) -h(y-z)\ < C + ez k (A.22) 

so that, altogether, 

H{y) - H(y - z) > -G + ezy^ 1 = -C + eu . (A.23) 
for z = tty 1- * . It is also easy to see that for each u > 0, 

lim [#(y) ~H{y~uy 1 ' k )] = aftli. (A.24) 

The estimate (|A.17)l follows then from the dominant convergence theorem applied to the integral on the 
right hand side of Eq. (|A. 18f) . 
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Appendix B 

We shall illustrate here the MFDT (14. 5p and its version (|4.11|) in the frequency space on the simple 
example of the one-dimensional Langevin equation (|l,12p on a circle with the Hamiltonian H = and a 
constant force G. In this case, Eq. (|1.12|l has, of course, the explicit solution 



x t = x + Gt + a/2/3" 1 W(t) (B.l) 

for the standard Brownian motion W(t). Since x is the angular variable, the above process possesses an 
invariant probability measure with the constant density g — The corresponding current j = ^G 
is constant and so is the local mean velocity v — G. The transition probabilities of the process have the 
form of the shifted and periodized heat kernel 

Pt(x,dy) = J— t }^ e 4t <^ = 57 /L e d y 

v n— 00 n— — 00 

= 27^(|(^-y + Gi),e- t//3 ), (B.2) 
where $3(2,5) is the Jacobi theta function [21]. For the (real) observables O a (x) = O n e mx , one 

n 

obtains easily the equalities: 

Tl ab (t-s) = 6(t - s)J20nO'l n n e {mG - p ~ ln2Kt - 8) , 

n 

C ab (t~s) = ^0 7 a l Oi n e lnG(t - s) - /3 ~ ln2|t - s| , 

n 

B ab (t - s) = -iG 6(t -s)J2 0»0-» n ^G-p-^^t-s) _ (B 3) 

n 

The result for the Fourier transforms (|4.9[) and (|4.10|l follows immediately: 

TZ ab (u) = ]To°O b _ n i - n ' — , 

' p L n* — iu) — inG 



2/T 1 n 2 



/3" 2 n 4 + (w + nG) 



2 ' 



S" b H = -J20 a n b . n 2 inG . (B.4) 

z — ' a 1 n 2 — ilo — inG 

n 

The modified FDT (|4.5|l or (|4.11|l are clearly satisfied. Taking a — b, one may define the factors 

A ft - s) = - - - — - for i>s, A (lu) = : (B.5) 

0d a C aa (t-8) (3cuC aa (uj) 

that control the violation of the standard FDT [131 114] . In particular, for O a (x) equal to sin(na;) or 
cos(na;) with one obtains (replacing the superscript a by n in this case): 

X (t - s) = [1 + — tan (nG(t - .)) j , X ( w ) = ^ + ^ 2 + ^ . (B.6) 

Note that X™(0) = 1 = X n (oo) but that these factors are not necessarily positive. This is also true for 
the "effective temperatures" 

T: s {t-s) = /3- X X a {t-s)-\ t: s {lo) = P^X^lo)- 1 . (B.7) 

In particular, T" b (lu) is positive only in the region where ui 2 > (n 2 G 2 — f3~ 2 n 4 ) and in this region it 
decreases with lu 2 approaching for u 2 — > 00 the value /3 _1 . 

The above calculations show that the equilibrium relation (| 1.2) 1 between the response and correlation 
functions is strongly violated in the Langevin equation on a circle with a constant drift unless the drift 
vanishes. On the other hand, the drift may be removed altogether by passing into the frame moving with 
constant velocity p — G, see Eq. (|B.1[I . This is captured by our MFDT 1)1.30 . Indeed, the solutions of 
the advection equation 1 )1.4) 1 are 



O a {t,x) = ]T O n e n{x - Gt) (B.8) 

n= — 00 

so that the Lagrangian- frame response and correlation functions are 

n ab {t, s) = n ab {t - S )| G ^ , ci% s) = c ab {t - S )\ G = (b.9) 

and the MFDT (fOJ takes the form of the equilibrium FDT holding for G = 0. 
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